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Abstract. We have performed a detailed study of the 
broadband spectra of the luminous (<; 10 36 erg s" 1 ) glob- 
ular cluster X-ray sources using BeppoSAX. With the ex- 
ception of X 2127+119, located in NGC 7078, all the other 
spectra are well represented by a two component model 
consisting of a disk-blackbody and Comptonized emission. 
The measured low-energy absorptions are in good agree- 
ment with those predicted from optical measurements of 
the host globular clusters. This implies that there is little 
intrinsic X-ray absorption within the binaries themselves, 
and that the above spectral model provides a good repre- 
sentation of the low-energy continua. The sources can be 
divided into two groups. In the first group, composed of 
3 ultra-compact (orbital period <1 hr) sources, the disk- 
blackbody temperatures and inner-radii appear physically 
realistic and the Comptonization seed photons tempera- 
tures and radii of the emission areas are consistent with 
the disk temperatures and inner radii. For all the other 
sources, the disk-blackbody parameters appear not to be 
physically realistic and the Comptonization parameters 
are unrelated to those of the disk-blackbody emission. If 
this is a spectral signature of ultra-compact binaries, this 
implies that no other ultra-compact binaries are present 
among those studied here. It is unclear why this difference 
between the two types of binaries should exist. One pos- 
sibility may be related to the mass ratio, which is similar 
in the ultracompact systems and binaries containing black 
holes. In the latter systems the soft components are also 
well-fit with disk-blackbody models, which appear to have 
physically realistic parameters. 



Key words: Accretion, accretion disks 
- Globular clusters - X-rays: general 



Stars: neutron 



1. Introduction 

The nature of the X-ray sources in globular clusters has 
provoked great observational and theoretical interest. Sur- 
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veys of the dozen brightest (<;10 36 erg s" 1 ) globular clus- 
ter X-ray sources have shown that they share many simi- 
larities with the low-mass X-ray binaries (LMXBs) found 
in the rest of the Galaxy - for example the presence of X- 
ray bursts and a lack of coherent pulsations in their per- 



sistcnt emissio n (He rtz & Wood |1985| ; Parmar et al. |1989 



Verbunt et al. 1995| ) . The presence of Type I bursts indi 



cates that the luminous systems contain accreting neutron 
stars rather than black holes, and they are generally con- 
sidered to be similar systems to the galactic X-ray burst 
sources. Their formation is, however, about a factor ^200 
more efficient inside globular clusters than in the galaxy 
itself. 

LMXB X-ray spectra are generally modeled by com- 
bining one or more spectral components. At luminosities 
~10 37 erg s _1 these often consist of a modified blackbody 
and some form of power-law like component thought to re- 
sult from the Comptonization of cooler photons. At lower 
luminosities, the spectra appear to simplify and are often 
fit with a power-law with an exponential cut-off at high 
energies (White et al. |l988|) . 

In this paper we present the results of a systematic 
survey of the luminous globular cluster X-ray sources un- 
dertaken as part of the BeppoSAX Core Program. We pay 
particular attention to the spectral properties of the per- 
sistent emission. A long-term goal of such studies is to 
understand how the spectral shape may be used to derive 
information on physically interesting parameters such as 
the mass accretion rate, the neutron star mass, radius and 
spin-rate, the system inclination and the viscosity of the 
material in the accretion disk. For surveys of LMXBs, the 
study of those located in globular clusters provides a num- 
ber of advantanges: (1) the distances to the many globu- 
lar clusters have been independently determined through 
main-sequence modeling and other methods (see e.g., Car- 
retta et al. 2000). This allows more accurate estimates of 
luminosity than is typical for galactic LMXBs. (2) The 
mean cluster abundances have been determined, allow- 
ing any dependence of X-ray properties on abundance to 
be examined. (3) The cluster reddenings have been mea- 
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sured so providing an independent estimate of low-energy 
absorption. For a number of clusters the reddening is sig- 
nificantly less than typical for galactic LMXBs. 

Results of the BeppoSAX observations of Terzanl, 
Terzan2, NGC 64 41, and NG C 7078 are to b e fou nd in 
Guainazzi et al. (Il999| , |l998j ), Parmar et al. fll99S| ) and 
Sidoli et al 



(200C) 



respectively. In addition, NGC 6440, 
Terzan6, and Liller 1 were observed by BeppoSAX in 
other programs. Results are reported in in't Zand et al. 
( |1999| , |2000D and Masetti et al. ( |2000| ), respectively. For 
these sources data were extracted from the SSD Bep- 
poSAX archive and reprocessed using the SAX Data Anal- 
ysis System (SAXDAS). 

2. Globular Cluster Properties 

Table [l] includes informations on the globular clusters 
which are known to contain luminous X-ray sources and 
have been analysed here. Most of the cluster distances, 
metallicities and extinctions are taken from Djorgovski 
( |1993D - In the case of Terzanl more recent estimates of 
the distance and extinction obtained using HST are used 
(Ortolani et al. 1999). For the highly absorbed globular 
cluster Terzan6 the results of Barbuy et al. (1997) are 
used, while the metallicity of Liller 1 is taken from the 



infrared measurements of Frogal et al. (1995). 

The metallicities of the studied globular clusters 
range from strongly underabundant compared to solar 
(NGC 7078, [Fe/H] = -2.17), to around solar (Liller 1, 
[Fe/H] = +0.25). Such a wide range of metallicity should 
allow any dependence of the X-ray properties on this pa- 
rameter to be probed. However, we caution that the ma- 
terial being accreted does not necessarily have the same 
metallicity as the host globular cluster since the compan- 
ion star may have undergone a non-standard evolution 
during which its envelope composition was altered. 

The extinctions to the studied globular clusters range 
from low (A v = 0.1, NGC 1851) to strongly absorbed 
(A v ~ 10, Liller 1). The amount of extinction strongly af- 
fects the low-energy cut-off of the X-ray data and so the 
sensitivity to any low-energy spectral components. As an 
example, the X-ray source in NGC 1851 is visible using 
BeppoSAX down to 0.1 keV, while that in Liller 1 is only 
detected above 1 keV. 

A summary of all previous observations of the Globular 
Clusters luminous X-ray sources can be found in Sidoli et 
al. ( |2000h| ). 



3. Observations 

Results from the Low-Energy Concentrator Spectrometer 
(LECS; 0.1-10 keV; Parmar et al. ~ 



1997) 



Medium-Energy 
Concentrator Spectrometer (MECS; 1.8-10 keV; Boella 
et al. |l99"7| ), High Pressure Gas Scintillation P roport ional 
Counter (HPGSPC; 5-120 keV; Manzo et al. |l997| ) and 



Phoswich Detection System (PDS; 15-300 keV; Frontera 



et al. 1997) on-board BeppoSAX are presented. All these 
instruments are coaligned and collectively referred to as 
the Narrow Field Instruments (NFI). The MECS and the 
LECS are grazing incidence telescopes with imaging gas 
scintillation proportional counters in their focal planes. 
The non-imaging HPGSPC consists of a single unit with 
a collimator that was alternatively rocked on- and 180' 
off-source every 96 s during most of the observations. Dur- 
ing the observations in the year 2000 the collimator re- 
mained on-source. The PDS consists of four independent 
non-imaging units arranged in pairs each having a sepa- 
rate collimator. Each collimator was alternatively rocked 
on- and 210' off-source every 96 s during the observations. 

Table [l] lists the BeppoSAX observations used in this 
survey together with some basic information about each 
globular cluster (see Sect. ||). The observations of Terzan2 
and NGC 6712 were made prior to the failure of one of the 
MECS units. Data from the entire observations were accu- 
mulated, except for the following sources: (1) X 1730—335 
(the Rapid Burster, hereafter RB). BeppoSAX carried out 
a total of 7 observations, partly analysed by Masetti et al. 
( |2000D , which were extracted from the archive and ex- 
amined for the presence of bursts. Only two observations 
included sufficiently long quiescent intervals for this study 
(see Table §. (2) For the 5.7 hour dip source X 1746-371 
the dipping and bursting intervals discussed in Parmar et 
al. ( pJ99| ) were excluded. (3) The eclipse observed from 
the 1 2.4 hr eclipsing system X 1747-313 (In 't Zand et al. 
|2000D was excluded. (4) The bursts seen from X 0512-401, 
X 1724-304, and X 1745-203 were excluded. Table | also 
indicates the instruments that were used in the spec- 
tral analysis of each source. In the case of the faintest 
source studied (Terzanl, 0.014 MECS counts s _1 ) only 
the MECS was used, while for the slightly brighter Liller 1 
observations the LECS and MECS were used. If available, 
all the NFI were used for the brightest sources. 

Since previous analyses used different software and cal- 
ibration releases, all relevant data were re-extracted. As 
usual, good data were selected from intervals when the 
elevation angle above the Earth's limb was >4° for the 
LECS and >5° for other instruments and when the in- 
strument configurations were nominal. This was done by 
using the SAXDAS 2.0.0 data analysis package. LECS and 
MECS data were extracted centered on the source posi- 
tions using radii of 8' and 4', respectively, corresponding 
to about 95% of the total source photons. In the case of the 
RB extraction radii of 4' and 2' were used in order to min- 
imize the contribution from the 32' distant bright X-ray 
source 4U 1728—337. Similarly, data from the non-imaging 
instruments were not included in the RB fits. Background 
subtraction for the imaging instruments was performed 
using standard files; this is not a critical procedure for 
the bright sources studied here (with the exception of 
Terzan 1). Background subtraction for the non- imaging in- 
struments was carried out using the offset pointing inter- 
vals. In the case of observations were the HPGSPC re- 
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Table 1. BeppoSAX observations of luminous globular cluster X-ray sources. T is the MECS exposure and R the 
energy range used for spectral fitting. LMHP indicate whether data from the LECS, MECS, HPGSPC, and PDS 
respectively, were included in the fits 



Globular 


X-ray 




Porb 




BeppoSAX Observation 


T 


R 


Inst. 


d 


Cluster 


A v 


Cluster 


Source 




(hr) 




Start 


Stop 


(ks) 


(keV) 




(kpc) 


[Fe/H] 




NGC 1851 


X0512- 


-401 


<0.85 


2000 


Feb 23 12:44 


Feb 25 07:44 


72.2 


0.1-25 


LM-P 


12.2 


-1.29 


0.1 


Terzan 2 


X 1724- 


-304 




1996 


Aug 17 04:29 


Aug 18 05:05 


37.1 


0.3-100 


LMHP 


10.0 


-0.25 


4.0 


Liller 1 


X1730- 


-335 




1998 


Feb 27 03:30 


Feb 27 19:22 


16.2 


1.0-10 


LM — 


7.9 


+0.25 


~10 










1998 


Mar 02 19:57 


Mar 03 12:35 


29.8 


1.0-10 


LM — 








Terzan 1 


X1732- 


-304 




1999 


Apr 03 12:52 


Apr 04 15:29 


47.1 


2.0-10 


-M — 


5.2 


-0.71 


6.8 


NGC 6440 


X 1745- 


-203 




1998 


Aug 26 02:11 


Aug 26 14:36 


21.6 


0.4-100 


LMHP 


7.0 


-0.34 


3.5 


NGC 6441 


X 1746- 


-370 


5.7 


1999 


Apr 04 18:07 


Apr 05 20:23 


49.9 


0.25-30 


LMHP 


10.7 


-0.53 


1.3 


Terzan 6 


X 1747- 


-313 


12.4 


1998 


Sep 06 06:16 


Sep 06 19:02 


28.9 


0.5-50 


LMHP 


7.0 


-0.61 


7.4 


NGC 6624 


X1820- 


-303 


0.2 


1997 


Oct 02 06:45 


Oct 02 19:90 


22.3 


0.3-30 


LMHP 


8.1 


-0.37 


0.9 










1998 


Apr 17 04:14 


Apr 18 03:12 


37.5 


0.3-30 


LMHP 
















1998 


Sep 23 12:44 


Sep 24 15:30 


48.3 


0.3-30 


LMHP 








NGC 6712 


X 1850- 


-087 


0.3 


1997 


Apr 24 11:47 


Apr 25 11:09 


42.0 


0.3-50 


LMHP 


6.8 


-1.01 


1.4 


NGC 7078 


X2127+119 


17.1 


1999 


Nov 16 01:15 


Nov 17 00:50 


35.5 


0.1-100 


LM-P 


10.5 


-2.17 


0.2 



mained coaligned, background subtraction was performed 
using data obtained when the instrument was looking at 
the dark Earth. 

The overall spectra extracted as described above were 
first investigated by simultaneously fitting data from as 
broad a band as possible. The LECS and MECS spectra 
were rebinned to oversamplc the full width half maximum 
of the energy resolution by a factor 3 and to have ad- 
ditionally a minimum of 20 counts per bin to allow use 
of x 2 statistics. The HPGSPC and PDS spectra were re- 
binned using the standard techniques in SAXDAS. Data 
were selected from within the energy ranges 0.1-8.0 keV 
(LECS), 1.8-10 keV (MECS), 8.0-20 keV (HPGSPC), and 
15-100 keV (PDS) whenever a sufficient number of counts 
were obtained (see Table |l|) . The photo-electric absorption 
cross sections of Morrison & McCammon (jl983 ) and the 
solar abundances of Anders & Grevesse (1989) are used 
throughout. Factors were included in the spectral fitting 
to allow for normalization uncertainties between the in- 
struments. These factors were constrained to be within 
their usual ranges during the fitting (Fiore et al. 1999| ). 
An additional uncertainty of 1% was added quadratically 
to the uncertainties to account for systematic uncertain- 
ties in instrumental responses. 

Since we wish to investigate the spectral properties 
of the sources, we first attempted to find a general spec- 
tral model that could be successfully applied to all the 
sources. Many luminous LMXB are fit with a two compo- 
nent model consisting of a blackbody, or blackbody-like 
component which might represent emission from an op- 
tically thick accretion disk or from the neutron star sur- 
face, together with a Comptonized component which may 
be interpreted as emission from a hot inner disk region or 
a boundary layer between the disk and the neutron star. 



For the soft component we tried both a single blackbody 
model and a disk-blackbody component, in the formula- 



tion of Mitsuda et al. ( |1984| ) and Makishima et al. ( pL986| ). 
The disk-blackbody model assumes that the gravitational 
energy released by the accreting material is locally dis- 
sipated into blackbody radiation, that the accretion flow 
is continuous throughout the disk, and that the effects 
of electron scattering are negligible. There are only two 
parameters in this model, ri n (cosi) ' 5 where r; n is the in- 
nermost radius of the disk, i is the inclination angle of the 
disk and kT;„ the blackbody effective temperature at r; n . 
Some limitations of this model are discussed in Merloni et 
al. ( p000| ). 

For the Comptonized component the xspec model 
COMPTT (Titarchuk 



1994; Hua 



Titarchuk 1995 



Titarchuk & Lyubarskij 199!:), which self-consistently cal- 
culates the spectrum produced by the Comptonization of 
soft photons in a hot plasma, was used. This model con- 
tains as free parameters the temperature of the Comp- 
tonizing electrons, kT e , the plasma optical depth with re- 
spect to electron scattering r p and the input temperature 
of the soft photon (Wien) distribution (kTo). A spherical 
geometry was assumed for the comptonizing region. 

Fit result to the broad-band spectra with the sin- 
gle blackbody plus Comptonization and disk-blackbody 
plus Comptonization are reported in Tables ^ and |3| re- 
spectively. A single temperature blackbody plus a cut- 
off power-law model was also tried, but with unaccept- 
able results in the majority of the sources. Since the disk- 
blackbody plus Comptonization emission (Table |^) always 
gives a better (or equally good) fit than single blackbody 
plus Comptonization emission (Table |^) , we will refer to 
it as the "standard" model. 
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Fig. 1. The BeppoSAX spectra of each of the globular cluster X-ray sources obtained using the standard model, 
in order of increasing 0.1-100 keV luminosity. The contributions of the two components (where evident) are shown 
separately. Except for NGC 6441 the COMPTT is the one extending to higher energies 
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Table 2. BeppoSAX spectral results using the COMPTT and single-blackbody model. For globular clusters with multiple 
observations, the results are given in the same order as in Table pL L is the 0.1-100 keV (assumed isotropic, in units of 
10 35 erg s _1 ) source luminosity using the distances given in Table]!] and assuming that the measured Nh is interstellar. 
The units of Nh are 10 22 atom cm~ 2 . fbb is the absorption corrected 0.1-100 keV ratio of blackbody to COMPTT fluxes. 
All uncertainties are quoted at 90% confidence 



Globular 


L 


N H 


kT 


kT e 




kTbb 


Rbb 


fbb 


X 2 /dof 








^jve v i 


l Kev ) 




[Ke v j 


\ Km ) 








NGC 1851 


43 


o.oo3±8:88I 


o.2o±8:8§ 


2.4±8 : ? 




o.56±8:8l 


6.7±?;f 





14 


153.2/97 


Terzan 2 


200 


0.83 ± 0.06 


i ii+ 015 


36±£ 10 


2-8±2!4 


o.66±8:8£ 


io.7=f=8:lo 





15 


197.1/178 


Liller 1 


8.7 


i.4±J:2 


<4.2 


2.5±g;? 


17=1=33 


<1.5 


<4.0 


<0 


14 


100.7/90 




2.5 


i.3±S:S 


0.43±g:g 


9.1±S 


7.3±Ii 










71.0/57 


Terzan 1 


0.059 


<0.66 


<0.48 


>2.2 


<20 










18.4/11 


NGC 6440 


66 


0.38 ± 0.06 


o.50±8:84 


9.6 ±0.1 


8.6±8:8o 


i.3±8:? 


o'.87±8;i 





05 


101.5/107 


NGC 6441 


140 


o.2o±g:gg 


0.51±8$ 


2.35±8!ol 


i7.2±8:?g 


o.27±8:S 


28±?1 





04 


117.6/106 


Terzan 6 


59 


1 Q4-0-01 


i.34±8:S 


16.0±?8 4 


2.o±8:?° 


0.65 ± 0.05 


6.7±1:S 





29 


144.5/119 


NGC 6624 


420 


o.io±8:o3 


0.38 ± 0.02 


2.97±8:8I 


13.3 ±0.2 


o.2i±8:8 5 


88 ±2 





05 


268.2/120 




500 


<0.015 


o.38±8:o 3 


o n + - 01 


13.3 ±0.1 


1 c;4- - 01 


1 Q4-0-07 





05 


232.5/117 




560 


o.25±8:8? 


o.o7±8:81 


O 7J.0.15 
' ^O.Ol 


15.8±S:g 


o.64±8: f 


14.4±8'7 





09 


588.4/134 


NGC 6712 


13 


o.i7±8:8! 


o.60±g : gg 


91±J8 


1 n , 0.05 


o.43±8:8? 


8.9±hi 





34 


181.1/154 


NGC 7078 


90 


o.2o±8:So 


o.o42±8:818 


37±i 


2.9±+^ 


1.1 ±0.05 


n q i_0.80 
°- >3=c 0.3Q 





24 


375.5/124 



Table 3. BeppoSAX spectral results using the standard COMPTT and disk-blackbody model. For globular clusters with 
multiple observations, the results are given in the same order as in Table [|. L is the 0.1-100 keV (assumed isotropic, 
in units of 10 35 erg s _1 ) source luminosity using the distances given in Table and assuming that the measured Nh 
is interstellar. The units of Nh are 10 22 atom cm~ 2 . f is the absorption corrected 0.1-100 keV ratio of disk-blackbody 
to COMPTT fluxes. 8 is the inclination angle of the disk. All uncertainties are quoted at 90% confidence 



Globular 


L 


N H 


kT 


kT c 




kT in 


ri„(cos6») - s 


f 


xVdof 


Cluster 






(keV) 


(keV) 




(keV) 


(km) 






NGC 1851 


47 


o.o26±8:8§8 


o.33±8:83 


O Q4-0.20 
^•°±0.50 


i4.o±8:g 


o.30±8:§I 


17±! 9 8 


0.14 


126.1/97 


Terzan 2 


220 


0.78 ± 0.05 


0.63 ± 0.02 


27±l° 


3.8 ±0.7 


3.3 ±0.3 


0.28 ±0.05 


0.12 


191.5/178 


Liller 1 


9.6 


2.8±?:l 


o.44±8:!l 


3.o±J:i 


6.08±8:i 


<0.50 


<37,000 


<0.05 


100.6/90 




2.5 


i.3±?:? 


o.43±8;l 


9.1±S 


7.3±I:J 








71.0/57 


Terzan 1 


0.059 


<0.66 


<0.48 


>2.2 


<20 








18.4/11 


NGC 6440 


71 


0.47 ±0.07 


0.50 ± 0.04 


8.7±S:i 


9.9±?;I 


l.98±8;?| 


0.47±8:Jg 


0.14 


84.8/107 


NGC 6441 


140 


0.26 ±0.02 


o.54±8:8l 


oc_i_225 
ZO±24.5 


o.40±Jio 


2.82 ±0.05 


o.92±8:8! 


10 


107.3/106 


Terzan 6 


67 


1.39 ±0.08 


o.66±8:?8 


8.3±li 


6.3±§;I 


2.65±8;o6 


0.60±8;8I 


2.16 


136.5/119 


NGC 6624 


480 


U.10±() 003 


o.574±8:881 


z - ' z 0.01 


1 c q_i_0.05 


o.57±8:8I 


22 ±0.13 


0.36 


152.1/120 




570 


o.i6±8:8o7 


o.46o±8:8? 


2.86±8:8! 


13.8 ±0.1 


o.37±8:8oi 


45±|° 


0.18 


182.2/117 




660 


o.27±8:8i 


o.o639±8:828 


2.77±8: i 


15.7±^g 


1.06±8^i 


4.9±?;I 


0.12 


627.7/134 


NGC 6712 


19 


o.39±8:8 3 


o.8i±8:i£ 


71±L 


i-7o±?:So 


o.63±8:8l 


4 -O±0.75 


0.76 


160.0/154 


NGC 7078 


65 


o.o6±8:8a 


o.o36±8:88l 5 


92±|;t 


2.6±8i 


1.92 ±0.04 


1.1 ±0.009 


2.0 


284.8/124 



4. Spectral Fit Results: the Standard Model 

The results of fitting the standard model to all of the 
extracted spectra are given in Table |[ All spectral uncer- 
tainties and upper-limits are given at 90% confidence. The 
fits to all the spectra are acceptable with the exception of 
NGC 7078 where the x 2 is 284.8 for 124 degrees of freedom 
(dof), and the second and third observation of the source 
in NGC 6624. Since a different model is a better represen- 
tation of the NGC 7078 X-ray spectrum, probably due to 



the ADC nature of this source (see Sidoli et al. 200C| ), we 
exclude it from the discussion, though the results of the 
fitting procedures have been included in Tables || and |[ 
for completeness. We note that an acceptable fit (a \ 2 of 
100.0 for 109 dof) can also be obtained without the disk- 
blackbody component for NGC 6440. In the case of Liller 1 
and Terzan 1 the sources are too faint to derive accurate 
spectral parameters. Moreover Liller 1 was not sampled 
at energies greater than 10 keV due to the reasons ex- 
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plained in Sect. 0. We include results from these sources 
in Tables || and g for completeness, but do not consider 
them further. Fig. [I] shows the photon spectra of the re- 
maining sources derived using the standard model. This 
figure clearly illustrates the importance of the broad-band 
coverage afforded by BeppoSAX in clearly separating the 
two components. Bolomctric luminosities have been cal- 
culated using the energy range 0.1-100 keV and assuming 
that the observed Nh is interstellar. The results in Sect 4.1 
indicate that for the majority of sources this assumption 
is valid. 



4.1. Low-energy absorption 

We first compared the fitted values of low-energy ab- 
sorption, Nh, with the values predicted using the re- 
lation between the optical extinction, A v , and absorp- 
tion of N H [cm- 2 /A v ] = 1.79 x 10 21 derived by Predehl & 



.00 



Schmitt ( |1995| ). If there is no additional absorption present 
in the X-ray binaries themselves, and their low-energy 
X-ray continua have been properly characterized, then a 
good agreement is expected between the values derived 
from the globular cluster A v measurements, and those 
from the LMXB X-ray spectra. Fig. || shows the mea- 
sured values of Nh taken from Table || plotted against the 
optically derived values of Nh ■ The straight line shows the 
locus of points with equal optical and X-ray derived Nh 
values. 

In general, there is good agreement between the two 
absorption measurements. In the case of X 1745— 203, lo- 
cated in NGC6440 the X-ray derived N H is -0.75 that 
expected (a difference of 1.5 x 10 21 atom cm -2 ). This may 
suggest the presence of an additional soft X-ray compo- 
nent that, not being included in the spectral model, leads 
to an underestimate of the absorbing column. Alterna- 
tively, the overall X-ray continuum model might be incor- 
rect. We note that also the blackbody plus COMPTT model 
does not provide a good estimate of the colum n density, 
being even much lower. In 't Zand et al. ( 1999 ) find that 
the spectrum of X 1745—203 may also be modeled using 
either a broken power-law and blackbody, or a high-energy 
cut-off power-law, or a bremsstrahlung and blackbody. In 
these 3 cases, the derived values of Nh are in the range 
6.7-9.0xl0 21 atom cm" 2 , similar to that predicted from 
the A v measurement of 6.2 x 10 21 atom cm" 2 . In the case 
of X 1850-087 located in NGC 6712 the X-ray derived N H , 
of 3.9 x 10 21 atom cm -2 , is significantly higher than that 
derived from the optical of 2.5 x 10 21 atom cm" 2 . This 
may suggest the presence of extra absorbing material in 
this system. However, when fitting this spectrum with a 
single blackbody plus Comptonization, a much lower col- 
umn density is found (1.7 x 10 21 atom cm" 2 ). 
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Fig. 2. The measured values of Nh (Table ||) plotted 
against the optically derived values of Nh- The straight 
line shows the location of points with equal optical and 
X-ray derived Nh values. Here and in all other figures, the 
light diamonds mark the X-ray binaries known to harbor 
sub-giant companions (P or f, > 10 hr), the triangles the 
ultra-compact binaries (P or b < 1 nr )> while the squares 
indicate all other sources. The meaning of the symbols is 
the same in all the figures 

4-2. Disk-blackbody component 

The soft emission from luminous globular cluster sources 
is succesfully fit with a multicolor disk-blackbody compo- 
nent. If we interprete this model "literally", its normal- 
ization translates into a lower limit for the innermost ra- 
dius r in of the accretion disk, Ri n =ri n (cosi) - 5 . Except for 
NGC 1851 and NGC 6624 (where the inclination cannot 
be constrained), the values observed are generally low, 
thus requiring a quite high inclination angle of the bi- 
nary systems, for rj n to be larger than the last stable orbit 
for an accretion disk around a neutron star (6GM/c 2 ~9 
(M/ M Q ) km). We note that still higher inclinations would 
be derived if the Merloni et al. (200C) correction to the 
disk-blackbody model were taken into account. For several 
globular cluster sources, a confirmation of their high incli- 
nation angles comes from X-ray and optical observations. 
In particular, the presence of strong sinusoidal-like modu- 
lations in the X-ray flux (as in the case of the ADC source 
NGC 7078), partial X-ray eclipses of the central source by 
the companion star (Terzan 6), X-ray dips (NGC 6441) 
are all indicative of inclination angles >60-70° (Frank et 
al. [19871 ) . 

The temperature kTj n at the inner disk radius is 0.5- 
3.5 keV, with the three ultra-compact binaries (NGC 1851, 
NGC 6712, NGC 6624; P orb < 1 hr) displaying signifi- 
cantly lower temperatures (<1 keV) compared to all other 
sources (1.5-3.5 keV). The variation of kT; n and inner ra- 
dius of the accretion disk is displayed in Fig. |3[ while the 
variation with the disk luminosity, measured in the 0.1- 
100 keV energy range, is shown in Fig. ^ and Fig. ||. 



Globular cluster X-ray sources 



7 



1 o.o r 



> 



, Ter6 6441 
6440 



1 .0 r 



^66?4 



1000 E 



X 

00 
X 



1 .0 



0.0 



00.0 



Rin (km) 



Fig. 3. Variation of the disk-blackbody temperature kTj n 
and the lower limit to the inner radius of the accretion 
disk 
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Fig. 4. Variation of the lower limit to the inner radius 
of the accretion disk and the luminosity of the disk- 
blackbody component, in the 0.1-100 keV band (in units 
erg s" 



In the disk-blackbody model, the bolometric luminos- 
ity, Lboh °f the accretion disk can be related to the temper- 
ature kTi n and the mass M of the compact object (see e.g., 
Eqs. 9 and 11 in Makishima et al. |2000| ) as L bo i oc M 2 -T in 4 . 
A summary of the parameters derived from the disk- 
blackbody component is shown in Fig. ||, where the glob- 
ular cluster sources appear to be divided into two groups: 
the ultra-compact binaries, with kTj n <,1 keV, the inferred 
masses of which are compatible with that of a neutron 
star, and sources with kT in <;1 keV, located on grids with 
much lower inferred masses for the compact object. Since 
all these X-ray binaries (except Terzan 6) are known to 
harbor a neutron star, due to their Type I X-ray burst- 
ing activity, either the multi-colour disk-blackbody is not 
an accurate model for the soft part of the spectrum, or 
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Fig. 5. Relation between the disk-blackbody luminosity 
measured in the 0.1-100 keV energy range (in units of 
10 35 erg s _1 ) and the temperature kT; n at the inner accre- 
tion disk radius. Dash-dotted lines represent loci of con- 
stant mass for the neutron star (in units of solar masses), 
while the solid lines show loci of constant r], that is the 
luminosity normalized to the Eddington luminosity (as- 
sumed to be 1.5 xlO 38 erg s _1 for a solar mass compact 
object) 



their luminosity has been underestimated, possibly due 
to a very high inclination angle of the binary systems. 
Another explanation may be that the parameters derived 
in the framework of the disk-blackbody model cannot be 
related directly to physical quantities (see Merloni et al. 

pooo| ). 



4-3. Comptonization component 



Thermal Comptonization (Titarchuk 1994) provides a 
good fit to the high energy region of the globular clus- 
ter sources spectra. The relations between the different 
parameters derived from the spectral analysis are shown 
in Figs. 6-11. The Comptonizing plasma has tempera- 
tures kT c ranging from few keV up to the extremely 
high value of ^92 keV in the case of the ADC source in 
NGC 7078. The seed photons in all sources have tempera- 
tures kTo<l keV while the Comptonization parameter y, 
defined as y=4kT c r 2 /m c 2 , is in the range ~l-7, except 
for the very low value of ^0.03 for the source in NGC 6441, 
that might be not in the right regime, since y^Cl and t<1. 
The trend of y with the luminosity is shown in Fig. ^, sug- 
gesting that at least for some sources a correlation exists. 
The dependence of the plasma optical depths r and of the 
electron temperatures kT c from the total luminosities are 
shown in Figs. ^ and ||. 

Following in't Zand et al. ( 1999| ), an equivalent spheri- 
cal radius Rw of the emission area of the seed photons can 

be defined as R w = 3 x 10 4 d A /%^/(fcT ) 2 km, where 
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Fig. 6. The Comptonization parameter y versus the to- Fig. 7. The optical depth r of the Comptonizing 
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d is the distance (in kpc), f CO mptt is the comptonized flux 
(in erg cm~ 2 s^ 1 , corrected for the interstellar absorp- 
tion, 0.1-100 keV band) and kTo is the seed tempera-, 
ture in keV. Rw ranges from few km up to ^1300 km 
(NGC 7078) and seems to correlate with the total lumi- " 
nosity (see Fig. as Rw oc L 5 . i 

This might suggest that the seed photons originate 
from the boundary layer (the radial extent and height of 
which become larger with increasing accretion rates; e.g. 
Popham & Sunyaev 2000) rather than from the accretion 
disk. The comparison between Rw and the lower limit 
Ri n to the inner radius of the accretion disk is shown in 
Fig. |TJ. 

In this respect, the comparison between the temper- 
atures kTo of the seed photons and the temperature 
kTj n is interesting (see Fig. |ll|): among the globular clus- 
ter sources, only the ultra-compact binaries (located in 
NGC 1851, NGC 6624 and NGC 6712) display tempera- 
tures kTo consistent with the temperature of the inner 
regions of the accretion disk. 



4-4- Connection with the globular clusters metallicity 

The galactic globular clusters hosting bright LMXBs are 
found to be both denser (except NGC 6712) and more 
metal-rich (except NGC 7078) (see Fig. |lj), indicative of 
the fact that both the high density stellar environment and 
the metallicity could play a role in the formation mecha- 
nism of LMXBs, thought to be via tidal captures or "t wo 
plus one" stellar encounters (see Bellazzini et al. 1995 ). 

Besides the formation mechanism, it is interesting to 
investigate if the metallicity of the accreting matter (as- 
sumed to be the same of the hosting globular cluster) plays 
a role in the LMXBs spectral properties. 



1 00.0 r 



0.0 r 



1 .0 T 



,644 1 [jfer2 



pa 

1851 



6624 
A 



0. 



10 



1000 



100 

Total Luminosity 

Fig. 8. The electron temperature kT e versus the total lu- 
minosity (0.1-100 keV, in units of 10 35 erg s" 1 ) 



In Figs. 13-19 several parameters derived from 
the spectral analysis are plotted versus the metallicity 
([Fe/H]) of the hosting globular cluster. 

The positive correlation between disk temperature 
kTj n and the metallicity displayed in Fig. [l3] is opposite 
to the results of Irwin & Bregman (1999), who studied 
in the ROSAT energy band the dependence of the X- 
ray properties of LMXBs on the metallicity of their envi- 
ronment, and found that the spectra of LMXBs become 
softer as the metallicity of their environment increases. 
However, we note that our correlation has a quite large 
probability of 11% to be obtained by chance; moreover, it 
might be spurious also for the reason that the three ultra- 
compact binaries, where the accreting matter should be 
essentially helium, are not expected to show a trend with 
the metallicity of the hosting globular cluster. Note that 
the same effect may be present when considering the ap- 
parent correlation between the total X-ray luminosity and 
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Fig. 9. The radius Rw of the emission area of the seed 
photons versus the total luminosity in the 0.1-100 keV 
band (in units of 10 35 erg s _1 ). Sources with larger lu- 
minosities (i.e. accretion rates) display larger radii of the 
seed photons regions 
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Fig. 10. Plot of the ratio between the lower limit R; n of 
the inner radius of the accetion disk and the radius R\v 
of the emission area of the seed photons, versus Rw 



the metallicity (see Fig. 15, probability of chance correla- 
tion of 20%). Moreover, no correlation is evident within 
the much more numerous X-ray sources residing in the 
M31 globular clusters (Verbunt et al. 1984), in the same 
range of metallicity. 

An apparent anticorrelation of the temperature of the 
seed photons kT with the metallicity is present, for 
sources in globular clusters more metallic than 0.1 solar 
(see Fig. |l7|). For the electron temperature kT e and the 
optical depth of the comptonizing plasma r there are no 
obvious correlations with the metallicity of the accreting 
matter. 
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Fig. 11. Comparison between the temperature of the seed 
photons, as estimated from the Comptonization model, 
and the temperature kTj n at the inner radius of the ac- 
cretion disk. The three ultra-compact binaries, indicated 
by the triangles, have all a temperature of the seed pho- 
tons consistent with originating from the accretion disk, 
while all other sources display a much higher inner disk 
temperature 
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Fig. 12. Globular clusters concentrations "c" 
(c=Log(cluster tidal radius)/Log(core radius)) ver- 
sus their metallicities ([Fe/H]). The 12 galactic globular 
clusters hosting bright LMXBs (squares) are located in 
the upper right region, being denser and more metal-rich 
than the other galactic globular clusters (marked by 
crosses) (see Bellazzini et al. |1995| ) 

4-.5. Iron line emission 

All the sources analysed here do not display significant 
iron line emission, except for the source in NGC 6624, 
where a very weak narrow line is detected. The addition 
of a narrow K Q iron line to the standard model, resulted 
in the upper limits to the EW reported in Table |]. 
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Fig. 13. Variation of the disk-blackbody temperature 
kTj n with the globular cluster metallicity ([Fe/H]). The 
accretion disk temperatures become higher as the metal- 
licity increases. However this correlation may be spurious, 
since the three ultra-compact binaries, where the accret- 
ing matter is essentially helium, are not expected to show 
a trend with the metallicity 
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Fig. 15. Total luminosity (0.1-100 keV, in units of 
10 35 erg s _1 ) versus globular clusters metallicity ([Fe/H]). 
A possible correlation exists, with higher luminosities as 
the metallicity increases. However the same caution as for 
the possible correlation of kT; n with [Fc/H] should be re- 
minded (see Fig. O) 
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5. Discussion 

We have performed a systematic survey of the properties 
of the luminous X-ray sources located in globular clus- 
ters. A remarkable property of this sample of LMXBs con- 
taining neutron stars is their well known distance and in- 
terstellar absorption. This, together with the broad-band 
capabilities of the BeppoSAX Narrow Field Instruments, 
allowed us to properly model the 0.1-100 keV spectra, 
to accurately determine the X-ray luminosity for these 
sources and to study their general spectral properties in 
a range of luminosity spanning more than two order of 
magnitudes. 

The spectral fitting gave satisfactory results with our 
spectral "standard" model, consisting of a disk-blackbody 



plus a Comptonizcd emission. Other two-component mod- 
els, such as a cut-off power-law plus a single-blackbody or 
a single-blackbody plus Comptonization emission, almost 
always resulted in poorer quality fits. 

However, although our standard model provides a good 
fit to the spectra of almost all the LMXBs studied, several 
problems arise with the physical interpretation of the re- 
sulting spectral parameters, especially when dealing with 
the disk-blackbody component. 

Indeed, the resulting inner disk radii are too small, 
even for extreme inclination angles of the accretion disks, 
and, on the other hand, the inner disk temperature are 
high compared with other systems in the galactic plane, 
for which kT; n ~1.5 is usually found (e.g. Tanaka & Lewin 



1995 ). For comparison, this same spectral model has been 
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Fig. 17. The temperature kTo of the seed photons in the 
Comptonization model, versus the globular cluster metal- 
licity ([Fc/H]). An anticorrelation exists for globular clus- 
ters more metallic than about 0.1 solar 
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Fig. 18. The plasma optical depth r versus the globular 
cluster metallicity ([Fe/H]) 

Table 4. Equivalent Width of a narrow iron line (energy 
fixed at 6.7 keV) added to the standard model. Upper 
limits are at the 95% confidence level 
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Fig. 19. The radius R\y of the region of the seed pho- 
tons versus the globular cluster metallicity ([Fe/H]). The 
correlation in the metallicity range 0.1-1 solar originates 
from the anticorrelation between kTo an d the metallicity 
in that range 



applied to the BeppoSAX data of two X-ray bursters, 
Ser X-l and GX 3+1, finding a r in (cosi) 5 of 6.8±0.8 km 
and 2.8± 1 2 km, respectively, and a disk temperature, 



kT in , of 1.46 and 1.95 keV (Oosterbroek et al. |2000J). 

The only globular cluster sources for which the multi- 
color disk-blackbody model seems to reasonably represent 
the likely physical parameters of the disk are the ultra- 
compact binaries, located in NGC 1851, NGC 6712 and 
NGC 6624. In these three sources the disk-blackbody pro- 
vides reasonable disk radii, and lower inner disk tempera- 
tures. Moreover, since in these three sources the temper- 
ature of the seed photons kTo an d the radius Rw are in 
agreement with kTj n and r; n respectively, it is tempting to 
identify the inner disk as the region for the seed photons 
of the Comptonization emission. However, we note that 
the finding of a lower disk temperature kTj n for the ultra- 
compact binaries could be partly due to the low interstel- 
lar absorption (as estimated from optical observations) to- 
wards these sources, with respect to other globular cluster 
sources (see Fig. po| ). 

For sources where temperatures kTo were significantly 
different from kT in (NGC 6440, NGC 6441, Terzan 2 and 
Terzan 6), we tried also to fit their spectra with the "stan- 
dard" model linking the temperatures kTo and kTj n to- 
gether. The main result of this fit with respect to those 
reported in Table |^ is the values of the disk-blackbody pa- 
rameters: the temperatures kTo=kT; n were in the range 
0.5-0.8 keV, the ratio f was less than 0.15 for all the 
three sources (remarkable, since previous analysis of the 
NGC 6441 spectrum resulted in a ratio f=10) and the 
disk-blackbody radii became 1-8.5 km. However, this al- 
ternative fit resulted in worse x 2 except for NGC 6440, for 
which an equally good x 2 was obtained. 
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Fig. 20. The inner disk temperature kTi„ versus the Nh 
estimated from the optical absorption (see Table 1) 

Another problem with the standard model arises from 
the relative contribution of the disk to the Comptonized 
component. Indeed, the ratio f of the disk-blackbody 0.1- 
100 keV luminosity to the Comptonized one, should be 



~1 (Sunyaev & Shakura 1986). But we find f>l for the 
sources NGC 7078 (the ADC source; f~2), Terzan 6 (f~2) 
up to the very high value of NGC 6441 (f=10). On the 
other hand, when a single-blackbody is used instead of 
the disk-blackbody, the ratio f^f, is always less than 1 (see 
Table ||). NGC 6441 is an intriguing source: if a blackbody 
is used instead of the disk-blackbody, the Comptonization 
dominates the emission and the high energy region of the 
spectrum, whereas when a disk-blackbody model is used, 
the opposite is valid. 

For the LMXBs where the interpretation of the soft 
component as contributed by the inner disk is reasonable, 
it is possible to estimate the strength of the neutron star 
magnetic field. The total luminosity L to t of an accreting 
neutron star can be expressed as the sum of the lumi- 
nosity liberated in the accretion disk Ldisk and the re- 
maining luminosity dissipated at the neutron star radius 
R ns , where L to t=GM ns M/R ns and L disk =GM ns M/2R m 
(e.g., Priedhorsky 1986 ). R m is the magnetospheric radius 

37 2/7 £ km, where B=B 8 10 8 G 
R 6 10 km, M ns =2 M 2 M is the 
37 ergs' 1 



R m = 20B^ /7 Rj o/T M^ /7 L 



10 37 L37 erg s 1 and the param- 



is the magnetic field, R n 
neutron star mass, L to t 
eter £ accounts for the disk geometry of the accretion and 
is £~0.5-l, (Ghosh & Lamb |1979[ |1992[ ). Thus, the mag- 
netic field strength of a neutron star can be estimated 
from the luminosity of the disk. 

In Fig. |2l] a grid with different neutron star magnetic 
fields, calculated assuming R6 = l, Ma=l and £=1, is su- 
perimposed on a plot of the ratio of disk to total luminos- 
ity, versus the total luminosity (0.1-100 keV). The r atio o f 
Ldisk to Ltot is a measure of R ns /2R m (Priedhorsky 1986|) : 
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Fig. 21. An estimate of the neutron star magnetic field 
strength depending on the disk and total luminosities in 
the 0.1-100 keV, in units of 10 37 erg s" 1 . Straight lines 
are loci of constant magnetic field (0.1, 0.5, 1, 2, 5 and 
10 x 10 s G, increasing from top to bottom), assuming a 
2 M Q neutron star. A lower neutron star mass results in a 
shift of the magnetic field grid towards the bottom of this 
plot. Only sources with f<l are reported 



Assuming a 2 Mq neutron star, this ratio depends only on 
the magnetic field strength and total luminosity. The ratio 
of Ldisk to Ltot reaches a maximum of 0.5 when R m is equal 
to R ns ; since the sources in NGC 6441, NGC 7078 and 
Terzan 6 have an unusually high disk-blackbody emission 
(as estimated from the disk-blackbody model used here), 
that led to ratios above the maximum allowed of 0.5, we 
will not consider them. 

The magnetic field estimated in this way for the source 
X 1820-30 in NGC 6624 (2.1x10 s G) is in agreement with 
the value derived from the saturation of the kHz QPO fre- 
quencies (Campana 2000] ). In the case of the RB, we find 
a magnetic field (~8x 10 8 G) consistent, within a factor of 
2, with that estimated by Masetti et al. ( 2000 ) using the 
propeller effect model. 
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6. Conclusions 

We performed broad-band spectroscopy of the luminous 
X-ray sources located in a number of galactic globular 
clusters with BeppoSAX. The results can be summarized 
as follows: 

— All the spectra can be acceptably fit by the same 
two-component model, comprising high energy Comp- 
tonization emission together with a disk-blackbody 
accounting for the soft part of the spectrum, except 
for the source in NGC 6441, where the disk-blackbody 
dominates the high energy region. 

— In the framework of this "standard" model and as 
a result of the spectral fitting, the sources can be 
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clearly divided into two groups, on the basis of their 
disk-blackbody parameters: (1) the sources located in 
NGC 1851, NGC 6712 and NGC 6624, and (2) all the 
other sources. Indeed, the disk-blackbody model, if 
literally interpreted, provides a physically reasonable 
interpretation for the first group of sources, while for 
the second group the disk-blackbody spectral param- 
eters are not physically acceptable, displaying high 
inner disk temperatures and extremely low innermost 
radii of the accretion disks, too low even for edge-on 
binary systems. There is no obvious explanation for 
a factor of ~3-10 lower disk temperatures in the 
ultra-compact binaries with respect to all the other 
systems; higher neutron star masses in ultra-compact 
binaries would contribute in reaching lower inner disk 
temperatures kTj n , but cannot explain, alone, a factor 
of ~3-10 lower disk temperatures with respect to 
other binary systems. Indeed, for systems with similar 
disk-blackbody luminosities, the ratio of the neutron 
star masses M ns i/M ns2 cx(kTi n 2/kTj n i) 2 , thus a factor 
of ~3-10 in inner temperatures translates in a factor 
of 10-100 in neutron star masses. 



except for the source in NGC 6624. This trend 
could be spurious, since in the three ultra-compact 
binaries the accreting matter might be only he- 
lium. On the other hand, if this is real, it suggests 
that matter accreting in the ultra-compact binaries 
NGC 1851 and NGC 6712 is probably not pure helium. 

The total luminosity L (0.1-100 keV) seems to in- 
crease with the metallicity, but again this trend could 
be spurious, since it is mainly based on the behavior 
of the ultra-compact binaries. 

There is no clear evidence for a correlation of the 
Comptonization spectral parameters with metallicity. 

We caution that an observational bias could be present 
due to the galactic absorption in the direction of the 
globular cluster sources: lower disk-temperatures kTj n 
have been derived for sources with low interstellar ab- 
sorption, and these happen to be the ultra-compact 
binaries. 



For the first group of sources, the Comptonization 
emission has seed photons temperatures kTo and radii 
Rw consistent with their inner disk temperatures 
kTj n and radii for the accretion disks ri n (cosi) , 
respectively. This could be interpreted with Comp- 
tonization sustained and fed by the inner regions 
of the accretion disks. Interestingly, this first group 
of sources is composed of the ultra-compact binary 
systems. If this result is indicative of a spectral sig- 
nature of the ultra-compact binaries, we suggest that 
the ultra-compact binary located in NGC 6652 will 
display the same spectral behavior and that, among 
the sources here studied for which both the orbital 
period and the optical counterpart are unknown, no 
other ultra-compact binary systems are present. 

In all other sources (group 2), the temperatures of the 
seed photons (kTo) ar e significantly lower than the 
temperatures kTi n . 

There is no strong evidence for a correlation of Comp- 
tonization spectral parameters with total luminosity 
L (0.1-100 keV), except for the radius Rw of the 
emission area of the seed photons, for which Rw oc 
L 05 . 



The search for correlations of the spectral parameters 
with the metallicity of the hosting globular clusters led to 
the following results: 

— The temperature kTj n of the inner radius increases 
as the metallicity of the hosting cluster increases, 
and the opposite is true for the inner disk radius, 



In a literal interpretation of the disk-blackbody model, 
the luminosity of the accretion disk, for the sources with 
a ratio f<l, is used to derive an estimate of the neutron 
star magnetic field. We note that, for the two sources for 
which an independent estimate of the magnetic field exists 
(the RB and X1820-303, located in NGC 6624), there is 
a good agreement with the values derived here. 

On the other hand, problems exist with a physical 
interpretation of the spectral parameters of the disk- 
blackbody component, since we obtained resonable val- 
ues for the radii of the inner disks only in the case of 
the ultra-compact binaries. From the similarities between 
radii and temperatures of the inner accretion disk and the 
parameters of the "seed" photons in the Comptonization 
emission, we propose that in the ultra-compact binaries 
in NGC 1851, NGC 6712 and NGC 6624 the Componiza- 
tion is fed by the inner regions of the accretion disk, while 
this is not the case for all other sources (in NGC 6440, 
NGC 6441, Terzan 2 and Terzan 6). It is unclear why the 
disk-blackbody model gives physically reasonable results 
only in the case of the ultra-compact binaries. A possibil- 
ity may be related to the mass ratio, q, which is similar 
in the ultracompact binaries and binary systems contain- 
ing black holes. Indeed, physically reasonable parameters 
are usually obtained from disk-blackbody fits to the soft 
emission from black hole X-ray binaries (see, e.g., Tanaka 
& Lewin|1995D. 
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